Supercurrent survival under Rosen-Zener quench of hard core bosons 
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We study the survival of super-currents in a system of impenetrable bosons subject to a quantum 
quench from its critical superfluid phase to an insulating phase. We show that the evolution of the 
current when the quench follows a Rosen-Zener profile is exactly solvable. This allows us to analyze 
a quench of arbitrary rate, from a sudden destruction of the superfluid to a slow opening of a gap. 
The decay and oscillations of the current are analytically derived, and studied numerically along 
with the momentum distribution after the quench. In the case of small supercurrent boosts f, we 
find that the current surviving at long times is proportional to u^. 



Progress in the field of experimental cold atom sys- 
tems allows a controlled and direct access to the non- 
equilibrium physics of quantum many body systems. 
This becomes particularly exciting when the system is in 
the vicinity of a phase transition and is driven through 
it as in the case of the superfluid(SF)-insulator tran- 
sition [1], or the magnetic ordering transition [2]. In 
such situations, complex physics is often exhibited; a 
prominent example is defect formation in the ordered 
phase due to a fast quench [3], which is qualitatively 
understood using the Kibblc-Zurek mechanism [4] origi- 
nally proposed in cosmology to describe domain forma- 
tion during cooling of the universe. Especially interesting 
are low-dimensional quantum systems under a "quantum 
quench" , with a Hamiltonian driven through a quantum- 
critical point. This can be achieved in quantum gases 
confined in highly anisotropic traps and in optical lat- 
tices [5, 6]. For instance, it was recently demonstrated 
that when a system is driven from an ordered to a disor- 
dered phase, the order-parameter correlations experience 
'revivar[l, 7-10]. Also, in the special case of driving 
an off-critical Id system into criticality, Calabrese and 
Cardy [11] showed that the long-time behavior of corre- 
lation functions can be obtained from the correlations in 
the final critical state. Despite these advances, as well as 
qualitative understanding of 'revival' phenomena and the 
Kibble-Zurek mechanism, analytical and exact results in 
this field are scarce. 

In this manuscript we focus on a system driven out of 
criticality by the variation of an external field (as opposed 
to Ref. [11], where a system is brought to criticality). In 
practice, this is a generic case, but elegant general results 
as in Ref. [11] for the opposing case are by and large 
absent. Special cases that have previously been stud- 
ied analytically are the behavior of a dipole model of a 
Mott insulator in an external electric field when the field 
is suddenly changed from a "no dipole" state to a maxi- 
mally polarized state [12], and the dynamics of traversing 
spin chains between two phases [7]. A related work, Ref. 
[8] , describes bosons driven abruptly from a Mott to su- 
perfluid phase and showed collective oscillations of the 
superfluid order parameter with period proportional to 



the gap in the initial Mott state. In [10], the evolution of 
hard-core bosons (HCBs) undergoing a quantum quench 
was studied numerically. 

A natural question arising in the SF-insulator transi- 
tion regards the fate of supercurrents in the system. In 
this manuscript we describe in detail the evolution of su- 
percurrents under quenching a 1-d system of HCBs, the 
so called Tonks-Girardeu gas [13, 14], from superfluid to 
insulator. Such systems may be formed by bosons at 
low temperatures and densities [6, 15]. We calculate the 
current as a function of time, while concentrating on the 
long-time current survival rate. Our analysis reveals fast, 
Bloch like, oscillations that are superimposed on a decay 
and a slower envelope function. The surviving current 
is found to be proportional to v^, where v is the initial 
supercurrent. Note that the decay of supercurrents was 
considered before mainly in Ref. [9] close to the Mott- 
superfluid transition. In addition, we present numerics 
for the evolution of the momentum distribution after a 
quench in the presence of a supercurrent. 

Our study of the supercurrent quench dynamics re- 
lies on a generalization of the Rosen-Zener problem [17] 
to the context of HCBs. The Rosen-Zener problem de- 
scribes a spin evolving in a time-dependent magnetic 
field with a particular profile, where an x-y magnetic 
field (analogous to a gap) is turned on while a z-field 
remains constant. In contrast, Landau-Zener like dynam- 
ics describe a spin in a constant x-y field, with Bz swept 
through zero; thus it is more useful for describing travers- 
ing the system through a quantum critical point [7] . The 
intcgrability of the Rosen-Zener evolution allows us to 
probe sudden quench dynamics as well as the response 
to a finite quenching time and so goes beyond previous 
treatments, which have dealt with an abrupt quench. 

The system we consider has the Hamiltonian 

H = -wJ2iblb,+i + h.c.) + V{t) Y,{-m\h (1) 

with b\ a boson creation operator at site i which obeys 
(6|)2 = 0, thus imposing the impenetrability of the 
HCBs. Throughout, we work in units such that the lat- 
tice spacing is 1. V{t) is the amplitude of an externally 
applied potential, which tunes the system from its super- 
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we note that k couples to fc + tt (or equivalently k — tt 
since k is the same up to multiples of 2tt) . Thus wc can 
write the Hamiltonian as: 
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FIG. 1: The Rosen-Zener switch-on profile V{t) for the super- 
lattice. 



fluid phase at F = 0, to an insulating phase at 7^ 0. 
When V{t) has the so-called Rosen-Zener shape, this evo- 
lution can be solved exactly. The Rosen-Zener profile is 
given by: 



Vit) 



Vo 



t < 0, 



cosh ^ ' 

Vo; t > 0, 



(2) 



where T is the turning-on. or quenching time. 

As we shall see, the HCB problem is equivalent to an 
infinite system of spins precessing in a time-dependent 
magnetic field. The evolution of a spin in the z direc- 
tion under = V{t) given in Eq. (2) (fig. (1)) was 
solved exactly by Rosen and Zener [17]. This exact non- 
equilibrium evolution allows us to establish the state of 
the system at t = 0, from which the system is evolved 
by the final Hamiltonian. We make use of this evolu- 
tion to analyze the fate of a supercurrent, introduced at 
t ^ —00 (when the system is in the gapless SF phase and 
V{—oo) ~ 0) by choosing our initial state as the ground 
state of the boosted Hamiltonian —w ^(e^'^&|&i_|_i -I-/1.C.). 

The first step in our analysis applies the Wigner- 
Jordan transformation to the HCB operators [18], 



(3) 



with aj being a fcrmion operator associated with site 
j. In the "fcrmion" picture, the presence of the super- 
current is expressed by choosing the "Fermi sea" of the 
a fermions to be shifted in momentum. That is, the k 
occupation of the fermions is 



ns{k) ~ Q{—kp + 1^ < k < kp + ly), 



(4) 



where kp is the Fermi momentum. In the following we 
will assume half filling, so kp = 7r/2. The current density 
is given by: 



J 



L ^ 
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where L is the total number of sites. 

To find the evolution of the current we first address 
the evolution of the fcrmion operators introduced in Eq. 
(3). By rewriting the Hamiltonian as: 



Hq^- Efc<|x| 2wcos(fc)4afe 

jj V(t) / t , t N 

Hd = 2^k<M[alak+n + al^^ak) 



(6) 



with: 



Hk{t) = 2wcoa{k)a,-Vit)a^ 



(8) 



acting in the {o-\'^\+Tr} mode space, and explicitly 
breaking the problem into a product of noninteracting 
spin-systems. Thus solving the evolution of the HCB 
system is equivalent to solving for the time evolution of 
an infinite series of time-dependent two-level systems. 

For a given fc, the Schrodingcr equation generated by 
Eq. (8) for a fermion operator i/j^ (t) = s{t)al. +p{t)al._^_^ 
can be rewritten as a second order differential equation 
of the form: 



T/2 



COsh(f )2 



S + {■iiwcos{k) 



-tanh(-))5 
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for S(t) = g2«"'cos(fc)tg^^^^ g^^^ ^Yic same equation for 
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~2iw cos{k)t 



p{t), with W 



-W. 



This equa- 



tion, via a change of variable z — i(tanh(^) -I- 1), can 
be recast into a hyper- geometric differential equation. 
The solution, satisfying the initial condition |S'(— oo)| = 
1, |P(— 00) = 0| (recall |fc| < ^), is given in terms of the 
hyper-geometric function: 



S{z) = 2-Fi(a, -a, c; z) 

p(z)^-»yi(T^ 2-XT"°°'''' 

2Fi{l + a,l~a,l + c;z) 



(10) 



where [24]: 



a 



VoT 1 2iwTcos{k) 

TT ' 2 TT 



(11) 



We now use the solution at t = (by setting z ^ 1/2 
in (10)) as a boundary condition for the dynamics under 
the final Hamiltonian. 

The evolution of the fermion operators at times t > 
is obtained by diagonalizing Eq. (8) at fixed V = Vq- 
This results in: 

aiit) ^ AkaiiO) + Bkal^jO) (12) 
4+.W=^4+.(0)-i3fe4(0) (13) 
where A, B are given explicitly by the relations: 



Ak ~ cos (Ekt) + i cos 6* sin (Ek) 
Bk = i sin 6 sin {Ekt) 



with — TT < k < TT, 
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4i«2 cos-^ k 



(14) 
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FIG. 2; Current evolution for an abrupt quench, T = 0. We 
set super-lattice strength Vo/w = 0.7 and supercurrent boost 



and 9 is defined through 



tan 6* 



Vo 



2w cos k 



(17) 



Using Eq. (10) with z ^ 1/2 to find a[.(0), a[._|_^(0) and 
substituting in Eq. (5), we find the current as a function 
of time 



(18) 



2«^« Ili%^sm{ka){2\soAk-poBl\^ - l){nk - nk+^) 

where Uk is the fermion momentum (Fermi-Dirac) occu- 
pation: (a^Ofc/) = Skk'^k- The above equations (10)-(18) 
constitute a complete solution to the problem of the cur- 
rent evolution under a Rosen-Zener quench. Below we 
focus on some interesting special cases. 

In particular, we concentrate on an initial state with 
a super-current at half filling and with a short switching 
time, wTv <^ 1, which allows non-adiabatic transitions. 
In this case, to leading order near k 7r/2, we may take 
c ~ 1/2, r(c) » r(l/2) = V^, giving: 



p{0) 



'VaT\ 



cos (Ml) 



(19) 



Assuming small supercurrents {v ^ 1), the current evo- 
lution, Eq. (18), simplifies to: 



2awi' 



where: A4{x,y) 



M{{2t-T)Vo,KvVt) 



cos{x)C{y) — sin{x)S{y) 



y 



(20) 
(21) 



and S, C are the Fresnel functions S{z) 

C{z) = d^cos(f u^), and k = 

Fig. 2 shows an example of the current evolution for 
an abrupt quench. Generically, the evolution is charac- 
terized by fast oscillations superimposed on a slow en- 
velope which oscillates and decays. The time between 



successive maxima of the slow envelope can be obtained 
by looking at the extrema of the Fresnel integrals in Eq. 
(21). From dzC{z) = cos(z^7r/2), we see that the extrema 
are at z = ^/2n, with n integer. Thus the condition 
\/2n = \/tvK yields the envelope period t^: 



(22) 



The fast current oscillations are a consequence of Bloch- 
like oscillations in each of the 2-state systems (k-state 
fermions and their fc — tt backscattered partners) de- 
scribed by Eq. (8). These arise from the cos(2t — T)Vb 
terms in Eq. (20), and thus have period tb = The 
Bloch oscillations decohere over time due to the spread 
in frequencies of the 2-level systems, given in Eq. (16). 
Indeed, we note that the current in Eq. (20) eventually 
decays, and no current is left in the system, to lowest 
order in v. 

Nevertheless, as the system is integrable and not cou- 
pled to a heat bath, some current survives at higher or- 
der in V. By computing the current averaged over long 
times, (j) = limr^oo ^ /q (j)dt, we show that the sur- 
viving current is proportional to for small supercur- 
rents. From Eq. (18), we find that the averaged current 
is given exactly by: 



_ Iw sin(i/) cos(VoT') 



1 - 



Vb arctan( 



^0 



2w sin(i^) 



which, to leading order in the supercurrent, is: 



■cos(\/or). 



(23) 



(24) 



which is of order , and thus is absent from the treat- 
ment leading to Eq. (20). 

To complement the current evolution analysis, we next 
analyze the momentum distribution of HCBs undergoing 
the SF- insulator quench while carrying a supercurrent. 
Although the analysis of the current evolution in the 
HCB system is the same for free fermions and a Tonks- 
Girardeu gas, the momentum distribution of the two sys- 
tems is quite different. While the presence of a supercur- 
rent is described by a shift of the Fermi step-function of 
the free Jordan- Wigner fermions, as in Eq. (4), in terms 
of the bosons, the supercurrent leads to a peak in the 
boson momentum occupation n/j = e'^'^'' {bj {t)bo{t)) at 
the boost value, i.e. Uk oc |/e — I'l"^/^ [14]. Upon quench- 
ing the system away from the critical superfluid, momen- 
tum modes at k undergo Bloch oscillations with modes 
at k + TT, resulting in oscillations of the current which 
eventually decays to an average, non-zero value. Unfor- 
tunately, an analytic description of the momentum distri- 
bution evolution is a much harder task, requiring detailed 
analysis of the determinants arising in the boson corre- 
lation functions. Even at equilibrium the analysis is not 
trivial, and relies on the translational invariance of the 
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FIG. 3: Uk as a function of k and time for a boost of = 
it/5, and super-lattice Vb = 4:W. The initial peak is at = 
V, and the peak transports back and forth to — 47r/5 before 
disappearing due to decoherence. Results obtained using an 
exact diagonalization study of chains 150 sites long. 




FIG. 4: Supercurrent surviving through a sudden quench. 
The symbols are the result of exact diagonalization of a chain 
350 sites long, (a) Current survival over the full range of 
boosts. From top to bottom, Vo/w = 0.1, 0.25, 0.5, 0.7 and the 
solid lines are Eq.(23). (b) The u"^ behavior of current survival 
at small currents. From top to bottom, Vo/w = 2,4,8. The 
solid lines are Eq.(24) 

system which allows the application of mathematical ma- 
chinery such as Szego limit theorems. Our case is much 
less accessible, due to the incoherent mixing between the 
different k modes. 

Using exact-diagonalization techniques (see e.g. [10]) 
to monitor the time evolution of the boson momentum 
distribution and the current, we investigated lattice sizes 
of up to 350 sites. The momentum distribution indeed re- 
flects the supercurrent Bloch oscillations, as can be seen 
in the bosonic nk plotted in Fig. 3. The period of these 
oscillations of the peak in between and h' — n agrees 
very well with the analytical result We also con- 
firmed numerically the main results of this manuscript, 
i.e., the current survival after a quench, Eqs. (23) and 
(24). Our results for the long-time averaged current (in 
a chain of 350 sites) vs. supercurrent boost v and super- 
lattice strength Vo/w for an abrupt quench are summa- 
rized in Fig. 4. Perfect agreement is obtained between 
the analytical and numerical results. 

In this paper we studied the behavior of hard-core 
bosons under a Rosen-Zener quench in the presence of 
a supercurrent. We described the evolution of the su- 



percurrent, and its long-time survival fraction, as well 
as the corresponding momentum distribution evolution. 
Perhaps the most readily accessible result is the V'^ law': 
starting with supercurrent v, the surviving current in the 
insulating phase is ex ly^ for i/ <g; 1. By using the Wigner- 
Jordan transformation, we essentially mapped the HCB 
gas to a Fermi-system with back-scattering at the Luther- 
Emery line [22]. In light of our results, it is particularly 
interesting to ask what happens when we consider a Lut- 
tinger liquid (describing either fermions or bosons): is 
the v^-la.w universal, or docs it depend on the Luttinger 
parameter, g < 27 This question, as well as a test of 
our predictions and the level of their universality can 
be taken on experimentally. This could be done, for in- 
stance, by probing cold-atoms in a ring-shaped trap, with 
a superimposed optical lattice. Such a geometry was re- 
cently discussed in Ref. [23]. If the additional optical 
lattice is made to rotate as it turns on, a supercurrent 
will exist in the rotating frame. The insulating phase 
can then be accessed by introducing a corotating super- 
lattice. Such an experimental setup will also be able to 
probe nonequilibrium quenches well beyond the regimes 
which we considered here analytically. 
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